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Method, apparatus and program for determining growth of 
trees 

BACKGROUND OF THE INVENTION 

5 1. Field of the Invention 

This invention relates to a solution for determining the growth of 
trees based on measurement data obtained at two different moments of time. 
Such measurement data can be obtained for instance by a laser scanner 
located above the trees. 

10 2. Description of the Prior Art 

Previous solutions for measuring the growth of trees are based on 
field measurements from the ground. Measurement equipment is thus 
transported to the location of a specific tree at two different occasions and the 
diameter or the height of the tree is measured from the ground. The obtained 

15 field measurement results are then compared to each other in order to 
determine the growth of the tree in question. 

The above described prior art solution has the drawback that it is 
very slow. The number of trees which in practice can be measured is relatively 
small due to the time needed and costs involved in carrying out 

20 measurements. 

Previously, for instance from WO-publication 01/31290, a solution 
is also known for obtaining a three-dimensional measurement result of stand 
attributes of trees located at a predetermined area. In this known solution, an 
aircraft, such as an airplane or a helicopter, is provided with measuring 

25 equipment which includes a laser scanner. Processing of the information 
obtained by the measuring equipment makes it possible to determine the size 
of the trees. This solution makes it possible to measure the size of a large 
number of trees during a relatively short period of time. So far, however, 
nobody has been able to disclose a solution for growth determination of trees 

30 where this known measurement solution is utilized. A significant problem has 
been the inaccuracy involved in the measurement. Present state of art 
methods do not utilize multi-temporal datasets in laser scanning. Even the use 
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of pixel-based change detection methods and multi-temporal laser data sets 
will result in inaccurate growth estimates. 

SUMMARY OF THE INVENTION 

-An object of the present invention is to solve the above mentioned 
5 drawback and to provide a solution which makes it possible to reliably 
determine the growth of trees significantly faster and in a more cost-effective 
way than in prior art solutions. 

Another object of the present invention is to provide a solution 
which makes it possible to process three-dimensional measurement data 

10 collected from an area at two different occasions such that measurement 
results can be utilized in order to calculate the growth of trees in this area. 

These and other objects of the invention are achieved with the 
method of independent claim 1, computer program of independent claim 6, 
and apparatus of independent claim 8. 

15 In this application the phrase 'tree location' refers to, for instance, a 

location where a tree, the crown of a tree, or parts of the crown of a tree are 
located. In case several trees are located very close to each other, the location 
of such a group of trees is interpreted as a tree location. In that case the 
determined growth at this tree location reflects the growth of the tree group. In 

20 case of such a group of trees, also the tree-to-tree matching is in this case 
carried out by matching the groups of trees. 

The invention is based on the idea that the measurement data 
collected at different occasions can be compared to each other when tree-to- 
tree matching is utilized. Tree-to-tree matching is carried out by, at first, 

25 identifying the tree locations separately from the first and second 
measurement data. Next, the locations, which according to both the first and 
the second measurement data are tree locations, are determined. The growth 
calculations are carried out only for such locations which are tree locations 
according to both the first and the second measurement data. This makes it 

30 possible to avoid errors resulting from positioning inaccuracy or from the Jact 
, that some of the trees in the area might have been cut or might have fallen 
between the moments of time at which the measurement data has been 
collected. 
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The most significant advantage with the present invention is 
considerable faster growth determination which make it possible to increase 
the number of trees whose growth is determined. 

In a preferred embodiment of the present invention, the average 
5 growth of trees is calculated by comparing the growth at a plurality of tree 
locations with at least one predetermined threshold value in order to identify 
tree locations where the growth is such that an error can be suspected, and 
calculating said average growth without taking into account the growth at said 
identified tree locations. This preferred embodiment makes it possible to filter 
10 out such tree locations whose growth is such that an error can be suspected, 
and thus to calculate the average growth by taking into account only reliable 
values. This improves the accuracy of the average calculations. 

Preferred embodiments of the method, computer program and 
apparatus of the invention are disclosed in dependent claims 2 - 5, 7 and 9. 

1 5 BRIEF DESCRIPTION OF THE DRAWINGS 

In the following, the present invention will be described in closer 
detail by way of example and with reference to the attached drawings, in which 
Figure 1 is illustrates collection of three-dimensional data, 
Figure 2 is a flow chart of a first preferred embodiment of the 

20 invention, 

Figure 3 is a flow chart of a second preferred embodiment of the 

invention, 

Figure 4 is a flow chart of a third preferred embodiment of the 
invention, and 

25 Figure 5 is a block diagram of an apparatus according to the 

present invention. 

DESCRIPTION OF PREFERRED EMBODIMENTS 

Figure 1 illustrates an example of a solution for collecting three- 
dimensional data to be used in growth calculations. 
30 In Figure 1, three-dimensional measurement data is collected by 

utilizing an aircraft 1 which flies a predetermined route over the area where the 
growth determinations are to be carried out. The aircraft 1 is provided with 
measuring equipment including a laser scanner 2. The laser scanner 
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comprises a laser gun for producing laser pulses and a detection unit, which 
records a received signal and determines the distance to a target. The time 
needed for the produced pulses to return to the laser scanner from the target 
is used for calculating the distance to the target. 
5 The laser scanner 2 can be for instance a TOPOSYS-1 -or 

TOPOSYS-2 laser scanner available from Topographische Systemdaten 
GmbH Wilhelm-Hauff-Strasse41. D - 88214 Ravensburg, Germany. Presently 
the measurements can be carried out for instance from an altitude in the range 
of 100 m to 1 ,5 km. If the altitude is for instance 800 m, the width of the area 

10 measured during one flight can be for instance 250 m. The measurements 
should be carried out such that laser beams hit the target substantially 
vertically. Tests have shown that incidence angles of more than 10° off-nadir 
result in a significant increase of shadowed areas. 

The orientation and position of the laser scanner 2 is typically 

15 defined with an inertia navigation system and with GPS (Global Positioning 
System) measurements. The inertia navigation system measures either 
orientation alone or both the orientation and position by using at least one 
inertia sensor. The GPS measurements are typically conducted using one 
GPS receiver 5 in the aircraft and another GPS receiver 6 on the ground as a 

20 reference station, usually within a 30 km range from the area where the 
measurements are conducted. 

The result of the measurements conducted with the solution of Fig. 
1 is three-dimensional measurement data, in other words a group of X, Y, Z- 
coordinates covering the entire measured area. This group of coordinates 

25 describes the shape of the area as seen from above. The number of 
measurement results obtained from an area naturally depends on the 
properties of the measuring equipment used. Practical tests have, however, 
shown that a density of about 3-10 measurement results per m^ is sufficient 
in order to identify individual trees from the measurement data. 

30 Figure 2 is a flow chart of a first preferred embodiment of the 

invention. The method illustrated in Figure 2 can be used to process three- 
dimensional measurement data obtained from the same area at two different 
moments of time. The collecting of data and the processing of the 
measurement data should be made exactly in the same way in order to obtain 

35 a result which is as reliable as possible. Thus, for instance, it is advantageous 
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if the aircraft could fly over the measured area by using the same flight paths 
both times. The algorithms used for processing the collected measurement 
data should also be exactly the same both times. 

The collected first and second measurement data is processed in 
5 blocks A and B. The processing is the same in both blocks, thus both the first 
and the second measurement data is in turns processed as will be explained 
in the following. 

A Tree Height Model (THM) is created as the difference between a 
Digital Surface Model (DSM) representing the top of the tree crowns and the 

10 Digital Elevation Model (DEM) representing the ground. This starts by 
generating the Digital Elevation Model (DEM) from the collected measurement 
data, for instance as described in WO-publication 01/31290. In this solution, 
the area covered by the measurement data is divided into pixels of suitable 
size, for instance 50 cm x 50 cm. The Digital Surface Model (DSM) of the 

15 crowns is obtained by taking the highest value of all hits with each pixel (50 
cm) and interpolating the surface model for missing pixels. In order to 
generate the Digital Elevation Model (DEM) from laser scanner data, the 
points that reflected from non-ground objects, such as trees and buildings, 
must be removed from the measurements. This can be accomplished by 

20 combining the high resolution Digital Elevation Model (DEM), which is 
generated by finding the minima of all the laser measurements belonging to 
each pixel area, with a low resolution one for identification and removal of non- 
ground points. This combination can then be repeated as often as necessary 
to generate a Digital Elevation Model (DEM) of the desired resolution. The 

25 procedure starts with a coarse resolution. Iteratively, the resolution is improved 
until the desired resolution is reached and a rough ground surface is created. 

The created surface is usually lower than a real surface because of 
the action of taking the minimum value for each DEM pixel, especially in hilly 
areas. Thus, a final refinement is performed to take return the original laser 

30 scanner data of the ground points and by comparing laser measurements with 
the corresponding values of the created surface and use the data in the final 
interpolation of DEM. 

The created Tree Height Model (THM) can be used for locating tree 
locations. The positions (X, Y coordinates) of the local Z maximum values (Z 

35 value describing the height of a tree at coordinates X, Y) represent the 
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locations of the tree locations. Preferably, the Tree Might Model is low-pass 
filtered before the local maximum search in order to reduce the number of 
detected tree locations. 

There are several methods to locate tree locations. Image 
5 processing techniques allow e.g. segmentation and crown delineation. yVith 
these knowledge you can calculate the center of the crown, the mass center of 
the crown and crown area. These techniques, also including morphological 
image analysis, are extensively reported in existing knowledge (books, 
scientific papers). In this new invention, all these techniques are applicable 

10 and do not change the basic novelty. All these techniques applied to 
distinguish tree locations are later on called as methods to locate the tree 
locations, for simplicity. The result is the tree locations according to the first 
measurement data (in block A) and the tree locations according to the second 
measurement data (in Block B). 

15 In block C, a search is carried out in order to determine such 

locations which are tree locations according to the processing of both the first 
and the second measurement result. If one of the processing results indicates 
that a specific location is a tree location but the other processing result 
indicates that the specific location is not a tree location, no growth calculation 

20 is carried out for that location. This tree-to-tree matching helps to ensure that 
growth calculations are carried out only for such locations which according to 
both processing results are tree locations, in practice the tree locations of the 
same tree or group of trees. 

The tree-to-tree matching can also be done with the help of crown 

25 areas. By overlaying the crown areas it is possible to calculate e.g. the 
distance between the mass centers of the crown area. This application 
assumes that the tree-to-tree matching by this mean is Included in the used 
term "determining the locations which are tree locations according to both said 
first and second processing results". 

30 The tree-to-tree matching carried out in block C eliminates from the 

growth calculations, for instance, tree locations were trees have been cut or 
were trees have fallen between the occasions when the two measurement 
results were obtained. 

Finally, in block D, the growth calculations are carried out for the 

35 locations which have been determined to be tree locations. The growth 
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calculation can be executed by calculating the height difference extracted from 
the two Tree Height Models (THM) for the given location. If the goal is to 
calculate the average growth in a given area, it can be carried out by 
subtracting the average of heights obtained for the two measurement results 
5 for all matched tree locations. 

The growth calculation can also be done in other ways. By using the 
image processing methods, e.g. segmentation to define the measured valuies 
or point clouds referring to each single tree (or group of trees), it is possible to 
compare the point clouds of the first and second measurement data. Since 

10 laser hits do not properly penetrate into crowns, the point cloud represents the 
. profile of each tree. By overlaying the profiles corresponding the first and 
second measurement, it is possible to reliable estimate the systematic shift of 
the profiles, which is due to horizontal or vertical growth of the trees. Also, 
depending on the laser point density on the trees and homogeneity of the 

15 stands, other statistical features (e.g. mean, median, histogram) can lead to 
improved accuracy in height determination. The growth calculations can thus, 
for instance, be done by calculating the difference of profiles, by calculating 
the vertical or horizontal difference of tree height models or by calculating the 
difference between the means of laser points for said first and second 

20 measurement data. 

Terrain elevations can be expected to be stable between the two 
occasions when the measurements were carried out. Therefore, the DEM 
created from two separate measurement results should result in the same 
elevation model. However, because Tree Height Model is the difference 

25 between the Digital Surface Model representing the top of the tree crowns and 
the Digital Elevation Model representing the ground, systematic errors will not 
show in the Tree Height Model, but random errors in DEM will propagate to 
Tree Height Model and therefore to the growth calculations. The effect of 
errors in DEM can be compensated in one of the measurement results as 

30 follows: 

- Systematic error of the two data sets can be determined for 
instance by measuring elevations of road surfaces. The measured difference 
can be compensated in one of the measurement results. 
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- The DEM value and difference between the two measurennent 
results can be calculated for each matched tree, and the growth can be 
compensated with the difference of the DEMs for each matched tree. 

Figure 3 is a flow chart of a second preferred embodiment of the 
5 invention. The solution of Figure 3 can be used in the tree-to-tree matching 
described in connection with block C of Figure 2. 

In block C1. the location, which according to one of the processing 
results (in this example the first processing result) is a tree location * is 
selected. 

10 In block C2, distance calculations are carried out in order to 

determine the minimum distance between the selected tree location and the 
location of the closest tree location according to the other processing result (in 
this example the second processing result). 

In block C3. the calculated distance DISTANCE is compared to a 

15 predetermined limit LIMIT. If the distance is smaller than the limit, it is 
determined in block C4 that both processing results indicate that the location 
in question is a tree location. Otherwise it is determined in block C5 that the 
location in question is not a tree location. 

The solution of Figure 3 has the advantage that the effect of small 

20 positioning errors regarding the tree locations can be eliminated. Setting the 
limit value LIMIT at 0,5 m, for instance, makes it possible to detect a tree 
match (a tree location according to both processing results) at a specific 
location even though there is a positioning error of 0,5 m in the measurement 
data. 

25 Figure 4 is a flow chart of a third preferred embodiment of the 

invention. This embodiment makes it possible to minimize possible errors in 
the measurement of the growth at single tree locations when the average 
growth calculations are carried out. This improves the accuracy of the average 
calculations. 

30 In block E a determined tree location is selected for comparison. 

Preferably, the process of Figure 4 is repeated until all tree locations have 
been systematically processed according to the indicated steps. 

In block F, the growth at the selected location is compared with at 
least one comparison value and in block G, a check is carried out to determine 

35 whether the comparison indicates an error. The object is to identify tree 
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locations where the determined growth is unrealistic. The number of 
comparison values and the value of these comparison values are selected 
with this in mind. An alternative is to use a first comparison value = 0, and to 
determine that an error exists for each tree location where the determined 

5 growth is equal to or less than zero. This is naturally a clear indication of an 
error, as it should not be possible that a tree becomes smaller by growth. A 
second comparison value can also be selected to Identify tree locations where 
the growth exceeds what is assumed to be realistic when taking into account 
the time period between the collection of the first and the second 

10 measurement data. For instance, it can be assumed that an error exists if the 
determined growth at a tree location is more than 1 m when the time period is 
one year. 

If it is determined in block G that the result of the comparison 
indicates an error, the tree location in question is not selected for average 

15 calculations but instead, a check to see whether all tree locations have already 
been processed is carried out in block I. 

Finally, when all tree locations have been processed, block J is 
entered. In block J, an average growth is calculated for the tree locations 
selected for average calculations in block H. 

20 Figure 5 is a block diagram of an apparatus according to the 

present invention. The apparatus 7 of Figure 5 can in practice consist of an 
ordinary Personal Computer (PC) programmed to carry out the growth 
calculations according to the method explained in connection with the previous 
Figures 2 to 4. 

25 The apparatus 7 has an input 11 for receiving first three- 

dimensional measurement data DATA1 and second three-dimensional 
measurement data DATA2. The data and the program used for processing the 
data is stored in memory 9. The processing means 8 carry out the processing 
and output the result of the growth calculations to a display 10. Alternatively or 

30 in addition to the display the result can also be outputted to some other 
computer peripherals, such as to a printer or written to a file in a disk. 

It is to be understood that the above description and the 
accompanying figures are only intended to illustrate the present invention. It 
will be obvious to those skilled in the art that the invention can be varied and 
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modified also in other ways without departing fronn the scope and spirit of the 
invention disclosed In the attached claims. 



